Non-adiabatic perturbations in decaying vacuum cosmology 
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We investigate a spatially flat Friedmann-Lemaitre-Robertson- Walker cosmology in which a de- 
caying vacuum term causes matter production at late times. Assuming a decay proportional to the 
Hubble rate, the ratio of the background energy densities of dark matter and dark energy changes 
with the cosmic scale factor as a~ 3 ^ 2 . The intrinsically non-adiabatic two-component perturbation 
dynamics of this model is reduced to a single second-order equation. Perturbations of the vacuum 
term are shown to be negligible on scales that are relevant for structure formation. On larger scales, 
dark-energy perturbations give a somewhat higher contribution but remain always smaller than the 
dark-matter perturbations. 



I. INTRODUCTION 



The origin of the accelerated expansion of the Universe, first described in is certainly the most challenging 
/"■s | problem of modern cosmology, receiving considerable attention today in both theoretical and observational areas 0). 

UThe simplest explanation is a cosmological constant A in Einstein's equations which remained the favored option until 
today and led to the ACDM model which also plays the role of a reference model for most studies in the field. Such 
, a constant is usually associated to the energy density of quantum vacuum fluctuations. Because of the cosmological 
Oh' constant problem in its different facets, including the coincidence problem, a great deal of work was devoted to 
q alternative approaches in which a similar dynamics as that of the ACDM model is reproduced with a time varying 
i-i ■ cosmological term, i.e., the cosmological constant is dynamiz ed. So far, the vacuum energy is not well determined by 
ryj I quantum field theories in curved space-time, given the divergences involved in its derivation. For the well-established 
renormalization method in the case of free conformal fields in an early de Sitter space-time 0, 0] the resulting energy 
density is proportional to i? 4 , where H is the Hubble parameter. But its value today is very small as compared to 
. the observed A. A correct order of magnitude is obtained for a dependence A ~ H 2 This circumstance gave 
J> ' rise to holographic dark-energy models (see, e.g. @, and references therein). Further semi-phenomenological ansatze 
C^l \ for a time-dependent cosmology term can be found in the literature Q- Among them there is the proposal that the 
QCD vacuum condensate associated to the chiral phase transition leads to a vacuum density proportional to H [8|, 
more precisely A = m 3 H, where m « 150MeV is the energy scale of the QCD phase transition. It is this approach 
that we are interested in in the present paper. The corresponding cosmological model is qualitatively similar to the 
ON standard one (the ACDM model) , with a radiation phase followed by a phase dominated by matter and, subsequently, 
by the cosmological term Q . The important point is that the time variation of the cosmological term is concomitant 
with a process of matter production in order to assure the conservation of the total energy. In this sense, the model 
may be seen as a particular case of interacting dark-energy models 0. The A oc H model has been tested from 
the observational viewpoint at both the background and perturbation levels. A preliminary joint analysis of the 
rcdshift-distancc relation for supernovas of type la, the baryonic acoustic oscillations and the position of the first 
peak in the anisotropy spectrum of the cosmic microwave background was performed in pd| . resulting in a present 
matter density parameter fiji/o ~ 0.36. At the perturbation level, the matter power spectrum has been calculated by 
assuming that the interacting cosmological term is strictly homogeneous fl2j . This spectrum reproduced the data for 
an even higher matter density parameter, fijv/o ~ 0.48. 

In the present paper we try to clarify the situation by reconsidering the previous analysis in a gauge-invariant setting. 
This allows us to understand the role of non-adiabatic perturbations in decaying vacuum cosmology. For simplicity 
we restrict ourselves to the dynamically most relevant components, namely dark matter and dark energy, i.e., we 
neglect both the baryon and the radiation components. We demonstrate that the entire two-component perturbation 
dynamics can be reduced to a single second-order perturbation equation for the density contrast of non-relativistic 
matter. We show also that at high redshifts the non-adiabaticity of the model is small, although non-vanishing, 
which allows us to use adiabatic initial conditions as a good approximation. In particular, we address the problem 
of the relevance of perturbations of the dark energy component for structure formation. In many investigations it 
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is assumed from the outset that dark energy does not cluster on small scales. But strictly speaking, this has to be 
justified on a case-by-case basis for all dynamical dark energy models. As it was argued by Park et al. [Hj|, neglecting 
the perturbations of the dark energy component may lead to inconsistencies and unreliable conclusions concerning the 
interpretation of observational data. For the present model of a decaying cosmologies! term, however, we shall show 
explicitly that the dark-energy fluctuations are indeed smaller than the dark-matter fluctuations by several orders 
of magnitude on scales that are relevant for structure formation. This justifies a posteriori the already mentioned 
previous analysis in which fluctuations of the dark energy component were not taken into account (l2j and which was 
compatible with a value of Ojv/o « 0.48 for the matter density parameter. This value is considerably higher than the 
corresponding ACDM- value and also higher than the background concordance value obtained in However, an 
update of the background tests with the most recent surveys of type Ia-supernovae seems to admit matter density 
parameter up to this value, at least for some of the data sets [l4j]. On superhorizon scales the fluctuations of the dark 
energy contribute a larger fraction to the total energy density perturbation but still remain smaller than the matter 
fluctuations. 

The paper is organized as follows. In section |n] we introduce the general expressions for our interacting two- 
component fluid, obtain the homogeneous and isotropic background dynamics of the model and perform a detailed 
analysis of the interaction term. In Scction lnll wc derive the central second-order equation for the matter perturbations 
in terms of which we also determine the perturbations of the dark-energy component. Furthermore, we discuss the 
non-adiabatic nature of the fluctuations and show that the non-adiabaticity is negligible at early times. Section IIVI 
presents the results of a numerical analysis which quantifies the role of perturbations of the dark-energy component. 
Finally, section [V] is devoted to conclusions and remarks. 



II. THE TWO-COMPONENT FLUID MODEL 



A. Basic properties 



We consider a two-component system with a total energy momentum tensor 

T ik = puiu k + ph ik , T% = 0, (1) 

where h ik = g$. + uiu k and gi k u l u k = — 1. The quantity u 1 denotes the total four- velocity of the cosmic substratum. 
Latin indices run from to 3. We assume a split of Ti k into a matter component (subindex M) and a dark energy 
component (subindex X), 

T ik =Tj% +Tf, (2) 

with (A = M, X) 

Tdk i k i t ik t ik ik , i k /o\ 

A = PAU A u A +p A h A , h A = g +u A u A . (3) 

Furthermore, we admit an interaction between the components: 

T^, k = Q\ Tf. k = -Q i . (4) 

Then, the separate energy-balance equations are 

- UMiT'u . k = PM,a.U a M + 6 A / (PM + Pm) = -U M aQ a (5) 

and 

- u X iT l x . k = px, a u a x + ©x (px +Px) = u Xa Q a ■ (6) 

In general, each component has its own four-velocity, with gi k u\u k A = — 1. The quantities <d A are defined as Q A = u a A . a . 
For the homogeneous and isotropic background we assume u a M = u x = u a . Likewise, we have the momentum balances 

Ki^m , k = (pm + Pm) u a M + PmM = h a Ml Q l (7) 

and 

h a Xl Tx ;fc = (Px +Px)u a x + VxM = -h a Xl Q\ (8) 
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where u A = u a A . b u h A . The source term Q % is split into parts proportional and perpendicular to the total four- velocity 
according to 



Q* = u l Q + Q i , 



(9) 



where Q = —UiQ 1 and Q l = h z a Q a , with UiQ 1 = 0. (Alternatively, one could have introduced a split with respect to 
the matter four-velocity. As we shall see later, for the model of interest here both options lead to identical results.) 

The contribution is supposed to describe some form of dark energy. In the simple case of an equation of state 
Px = —px-, where px is not necessarily constant, we have 



We are interested here in the case 



rpik ik 

1 x — -pxg 



px = 3 e , 



(10) 



(11) 



where = u? a is the expansion scalar and er is a constant. In the homogeneous and isotropic background one has 
O = 3H, where H is the Hubble rate. 



B. Background dynamics 



In the homogeneous and isotropic background we have 

3H 2 

and, assuming pm = from now on, 



Combining (fT2j) and (fT3|) we obtain 



nGp = 8nG (pm + px) = 87rG (pM + crH) 
H = -4:nG(p + p) = -AnGpM ■ 



(12) 
(13) 



H 



--H 2 +4mGaH . 



With 



3Hq = 8irGp , SImo = and a = (1 - Qmo) 

Po -Ho 



integration yields the Hubble rate @ 



H = H Q 



(14) 



(15) 



(16) 



where a subindex indicates the present value of the corresponding variable. The quantity a denotes the scale factor 
of the Robertson- Walker metric which was normalized to oq = 1. The energy densities are given by 0] 



Pm 
Po 



and 



^ = (1 - Omo) 
Po 



1 — f2jwo + ^A/oa 3 ^ 2 



(17) 



(18) 



respectively. For a = as well as for a << 1 we consistently recover the Einstein - de Sitter universe. For a >> 1, 
the solution tends to de Sitter space-time. The ratio of the energy densities is 



Pm 



PX 1 — ^MO 



MO -3/2 

a ' 



(19) 



It decays with a 3 / 2 , i.e., with a lesser rate than in the ACDM model for which the corresponding ratio decays as 
a -3 . But in both cases it approaches zero in the long-time limit. 
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The balances ([5]) and (|6|) take the forms 

p M + 3H PM = Q° (20) 

and 

Px = -Q° , (21) 
respectively. Since px is given as px — f @ = oH, the background source terms are 

Q° = u°Q = Q=px = ~Q = -a-H and Q a = . (22) 

Alternatively, 

u a Q a = -Q = -px = ^Q = (jH and Q a = . (23) 
Together with (fT5j) we find the balance equations 

p M + 3ffp M = AnGapM (24) 

and 

p x = -47rGcrpAf . (25) 
Differentiating the expressions (fT7|) and (|18p one realizes that they indeed are solutions of (|24p and (f2~5j) . 



C. The perturbed source term 

Denoting first-order perturbations by a hat symbol and recalling that for the background u M = u x = u a is valid, 
the perturbed time components of the four-velocities are 

u = u° = u° M = u° x = igoo • (26) 

According to the perfect-fluid structure of both the total energy-momentum tensor ([1]) and the energy-momentum 
tensors of the components in ([3]), and with it M = u x = u a in the background, we have first-order energy-density 
perturbations p = pu + Px , pressure perturbations p = pM + Px = Px and 

T a = T Mot + T Xa (P + P) ^ = pM^Ma + (PX + Px) U Xa ■ (27) 

Greek indices run from 1 to 3. For px = —px it follows 

px = -px =>■ p+P = Pm => u Ma = u a . (28) 

Since the component M is supposed to describe matter, it is clear from (|27[) that the perturbed matter velocity um<x 

coincides with the total velocity perturbation u a . With — u n up to first order, the energy balance in ([5} (correct 
up to first order) can be written as 

PM, a u a = —QpM - u a Q a . (29) 
On the other hand, the total energy balance is 

P , a u a = -e(p + P ) . (30) 

For the difference it follows that 

p- p M = {p- p M ) a u a = u a Q a . (31) 
Since, at least up to linear order, p — pu = Px, equation (f3~Tj) is equivalent (up to the first order) to 

Px = Px. a u a - u a Q a . (32) 
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In zeroth order we recover (|2 1 1) with (|22j) . The first-order equation is (cf. (|26|) ) 

Px + Pxu° = (u a Q a ) . (33) 



Notice that ([33| results from a combination of the total energy conservation and the matter energy balance. It has 
to be consistent with the dark energy balance ©. At first order, the latter becomes 

'p x + Pxu° = {u Xa Q a ) ■ (34) 

This means that 

(uxaQ a ) = (u a Q a j , (35) 
i.e., the projections of Q a along ux a and along u a coincide. Explicitly, 

(u a Q a j = (u a u a Q) = -Q = |e . (36) 
In a next step we consider the momentum balances. The total momentum conservation is described by 

K T % = (pm +Px+Px)u a + h ai Px ,i = . (37) 
With px = —px and px = — § B, we have 

h?T% = p M u a -h ai ^Q,i = . (38) 
Using u M = u n again, the momentum balance ([7]) for the matter component becomes 

K T M ;k = PMU a = h al Q t . (39) 

Consistency between ([38]) and (|39|) requires 

Q a = h ai Qi = ^h ai e,i . (40) 

Notice that we have only used the total momentum conservation and the matter momentum balance. The momentum 
balance ([8]) of the dark energy degenerates for the case px = —px- It does not describe any dynamics. 

With (|40[) . the source term in the momentum balance is explicitly known. Up to first order the matter momentum 
balance ([5^]) reduces to 

PMUa = Qa ■ (41) 

Explicitly, from (1401) the first-order source term becomes 

Q = o , Q a = | [e, Q + u a e\ . (42) 



D. Basic set of equations 



Before explicitly starting the perturbativc analysis, we summarize the basic set of equations. It consists of the 
energy balance 

Pm + &pm = -u a Q a = Q = ~Q , (43) 

the momentum balance 

PMUa = , (44) 

and the Raychaudhuri equation 

e + le 2 - u° a + 4nG ( p + 3p) = o . (45) 

o 

In the present case (j>m = and px = —px = — §©) the latter reduces to 

9 + ie 2 - u» + 4itG P m - ^-<?0 = • (46) 

Notice that the source terms in the balances (|43]) and (|44"]) are given by derivatives (temporal and spatial, respectively) 
of the expansion scalar. In the background, the energy balance (|43|) reduces to ([20| with (|22|) . the momentum balance 
(|44p is zero identically and the Raychaudhuri equation (upon using Friedmann's equation) specifics to (|14[) . The 
zeroth-order solutions of the system (|4"3")l . (|44[) and (|4^|) arc given by (fTT|) . (|T5)) and (JTHJ) , respectively. 
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III. PERTURBATION DYNAMICS 



A. General relations 



The general line element for scalar perturbations is 

ds 2 = - (1 + 20) dt 2 + 2a 2 F a dtdx a + a 2 [(1 - 2ip) 5 af} + 2E >a p] dx a dx p 
For the perturbed spatial components of the four- velocity one has 



a 2 rt L + a 2 F^ = Up = v tfJj , 



(47) 



(48) 



which defines the velocity perturbation v. A choice v = corresponds to the comoving gauge. Taking into account 
the definition pm = PM.a u<1 , the balance (|43[) becomes in first order 



S M + AmGa5 M - <p (-3H + inGa) + 6 = — , 

Pm 



(49) 



where we have introduced the first-order fractional perturbation 5 m = ~"- With the definition u a = u a] bU which 
leads to (u a ) = u a + 4>,at the momentum balance (|44[) becomes in first order 



o 



V + 



e + Qv 



3pM ^ ' ,a 3pM 

To calculate the source term Q in (|36[) we have the perturbed Raychaudhuri equation. With 



1 



ipMCL 1 



A6 + 6Aij 



where A is the three-dimensional Laplacian, one obtains 



e = _ 06H r 

3 3p a/ a 



CT / a A A a \ c- 87rG ~ 

A6 + 6Aw - A-kGpmSm + -^—aQ 



It follows that the source term -2- in (|4"9"j) becomes 



Q_ 

PM 



8ttG 2 \ 9 

a B 

3 3 



AB + BAw - 4ttG<5 



>M 



It is now convenient to describe the dynamics in terms of the gauge-invariant quantities 

B c = B + Sv , S C M = 5 M + ^% and 5 C X = 5 X + ^-w , 

Pm Px 



(50) 



(51) 



(52) 



(53) 



(54) 



which characterize the perturbations of the expansion scalar, the matter density and the dark energy density, respec- 
tively, on comoving hypcrsurfaces. Then, equation (|4l?]) takes the form 



J M 



4ttGct5 



M 



3pM 



(3H - AirGa) 



Pm 



(55) 



Here, Q c = Q + Qv is the gauge-invariant perturbation of the source term Q with the explicit (momentum space) 
structure 



Pm 



8ttG „„ a k 2 \ 6 C 
a - 2H 



3pm a 2 J p M 



AttGS 



M 



(56) 



where k 2 is the square of the comoving wave vector. Eq. (|55|) with (|56j) represents a combination of the energy and 
momentum balances of the matter fluid. Combining (|55p and (|56[) provides us with 



6 c M + ^a5 M + KQ c = 



S°M + 



8ttG 



M 



(57) 
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where 



K = 1 



aH 

3pM 



aH k 2 



4ttG 



3 H 3p M a 2 H 2 



(58) 



In the next step we differentiate (|57[) . eliminate C through the perturbed Raychaudhuri equation and C through 
(|57p and change to the scale factor a as independent variable. The result is the following closed second-order equation 
for 6 C M (a): 



(59) 



where the prime means a derivative with respect to a. The coefficients g(a) and f(a) are explicitly known and depend 
only on the background dynamics. They are given by 



9W 



and 



with 



f(a) 



2 K 



1 k 2 

-A- 



3 H 2 a 2 



3 3 
- — -B K 
2 2 



, L 1 1 , k 2 

[ K- 2 -2 B --, A ^H 2 



A = 



aH 

Pm 



i - n 



MO 



n M0 a- 3 / 2 ' 



(60) 



(61) 



(62) 



B = 



8ttG a 
~3~H 



i - n 



MO 



MO 



M0 a 



-3/2 



(63) 



K = 1 



3 A -6^-9 A ^ 



and 



1 1 9 - 



(3B - 4) 



A- 2 



a 2 H 2 



(64) 



(65) 



Equation (|59p for <5^ encodes the entire perturbation dynamics of the model, it is the central relation of the paper. 

with 8 C from (1571) . It is determined 

by S[ 



J M 

The perturbation 8\ = ^ of the dark energy can be obtained from p c x 
and its first derivative: 



3 w 



M 



1 



[aS 



BS 



Mi 



(66) 



The coefficients A and B do not depend on the wavenumber fc, they are of the order of one around the present time 
with a ~ 1. The quantity K, however, is scale-dependent. On scales well inside the horizon, equivalent to ^ ^ H 2 , 
one has \K\ ^ 1. This suggests that, on scales which are relevant for structure formation, the fluctuations 8 C X of the 
dark energy should be very small compared with the matter fluctuations S M . The quantitative analysis of section llVl 
below will confirm this behavior. It is expedient to note that Eq. ([5^1 is of the same structure as the corresponding 
perturbation equation for a unified viscous dark sector fluid (l5j . However, the coefficients g(a) and /(a) are different. 
In particular, the wave-number dependence is more complicated for the present model. 



B. Non-adiabatic perturbations 

To clarify the non-adiabatic character of the model it is useful to consider the gauge-invariant perturbations of the 
dark pressure, p x = px +Pxv and of the dark energy, p° x = px + pxv. Let us assume that in the corresponding rest 
frame p c x = c 2 p c x . Then, p x may generally be written as 

Px = c 2 s p x + (c 2 s -^-) vp x ■ (67) 
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In the general case, the non-adiabatic part of the perturbations of the X component is given by 

Px - —Px =\c a - 1 —\px- (68) 



Px \ Px 

For the present model c 2 s = = — 1 is valid, i.e., the X component by itself is adiabatic. For the total non-adiabatic 
perturbations we have 

- P - " PX - . . PX I - , * s , an , 

P~ ~P = Px P = Px + — [Pm + Px) ■ (69) 

P P P 

Introducing here (|67[) and rearranging we obtain 

_ -s P - PXPM f PM PX \ . / 2 , -, \ ~c rn n \ 

Pnad =p--p= ; — + {<% + 1) p x . (70 

P P \PM PX) 

Although the X component has no intrinsic non-adiabaticity for the present case, c 2 s = — 1 and px = —px, the 
two-component system as a whole is non-adiabatic. As it is obvious from ([TOl . non-adiabatic first-order perturbations 
can be characterized by the quantity 

^ - P2L + 0. (71) 



Since 



Pm Px 



Ph - 5 ° M and fk = (72) 



Pm —3H + AttGo px 127rGpA/ 

where we have used (|24j) and (|23|) . respectively, it follows with ([57]) that 

Pm Px _ s< m 2 Po 



Pm Px -3H[l-*ffi] 3 3H§ PM K 



7 [S M aH + H {l-n M0 )S c M }. (73) 



Notice that 



Pm Px = P M Px 
Pm Px Pm Px 



(74) 



i.e., the combination (|71|) is gauge- invariant, although the single terms by themselves are not. Equation (|73j) implies 
that the non-adiabatic contribution is completely determined by 5 M and its first derivative. Once 5 M (a) is known, the 
expression (|73[) is determined as well. The existence of non-adiabatic perturbations is characteristic for interacting 
two-component systems. In our case, the non-interacting limit corresponds to a = => px = =>■ Omo = lj 
equivalent to the one-component Einstein - de Sitter universe which, of course, has a purely adiabatic perturbation 
dynamics. 



C. Almost adiabatic initial conditions 



For sCl, the non-adiabaticity (|73]l can be calculated explicitly. This is relevant for specifying the initial conditions. 

It will turn out that, for the present model, almost adiabatic initial conditions are appropriate. At high redshifts, i.e., 
for a ^ 1, Eq. (|59|) reduces to 

S M + ^S M -^5 M = (a«l), (75) 

with solutions 

S C M = cia + c 2 a- 3/2 => S M = ci - ^c 2 a" 5/2 (a < 1) . (76) 
The asymptotic values for H, pm and K for a <C 1 are 

if ss tf ^Moa- 3/2 , Pm w Po^mo«" 3 , A' w 1 (a « 1) . (77) 
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Applying (JTBJ and dTTJ) in (73]) wc find, for a < 1 



Pa/ _ Px _ 
Pa/ Px 3-ff ^A/o 
2 a 3 



9 ff 2 S! 2 



3/2 

cia + C2a ' 

ci - |c 2 a- 5 / 2 ) a^^HoVtMo + H (1- fl M0 ) ( Cl a + c 2 a~ s / 2 



(78) 

The a-indcpendcnt C2-terms cancel each other exactly. The leading term for the growing mode is 

Pm Px » -"0" A/0 

It is interesting to note, that both terms on the left-hand side contribute with comparable amounts, although p c M 3> p c x . 
The non-adiabatic pressure perturbation (|70[) associated with the growing mode becomes 

5 Po^A/O (1 - f^A/o) , ^ , Qn , 

Vnad~-z c \ tb (a<l), 80) 

the corresponding fractional quantity is 

Pnad 5 1 - OmO 5/2 / 

— — - cia 0/ (a<l). (81 

p + p 6 S2a/ 

It is only for Qmo = 1? corresponding to the one-component Einstein-dc Sitter model, that the non-adiabatic con- 
tribution exactly vanishes. Although the perturbation 5 X in (|66[) is different from zero for Qmo = 1, it does not 
contribute to the total energy density perturbation since the background density of the X component vanishes (see 
relation (|84|) below). The ratio 

^ii-%V/ 2 (a«l) (82) 

\Pm\ a "A/0 

of the perturbations scales as the background ratio for these components (cf. Eqs. (fT7j) and (fT8|l). As expected, the 
non-adiabatic pressure perturbation is of the order of the dark-energy perturbation: 

I Pnad I 5 / „ -\ /oo\ 

(a < 1) . (83) 



IPX | 2 

It follows that at high redshifts it is not only the background contribution of the dark energy density which is small 
compared with the total background energy density, the dark-energy perturbations, together with the non-adiabatic 
pressure perturbations, are much smaller than the total energy-perturbations as well. 

On the other hand, Eq. ([75]) is also the correct one-component limit of Eq. ([5^1 for S7mo = lj i- e -i the Einstein- 
de Sitter model with a purely adiabatic perturbation dynamics. Consequently, the smallness of the non-adiabatic 
contribution at high redshift allows us to use approximately adiabatic initial conditions for k > 0. However, the mere 
two-component nature of our model requires that there is always a small, non- vanishing admixture of non-adiabaticity. 
This is consistent with the result in [16j , that purely adiabatic large-scale perturbations can never give rise to entropy 
perturbations. 

IV. DARK ENERGY PERTURBATIONS AND POWER SPECTRUM 

In the present study we have neglected baryons. This can be considered a reasonable approximation since baryons 
represent only 5% of the total amount of energy. Nevertheless one should keep in mind that the observed power 
spectrum reflects the distribution of luminous, baryonic matter. For the ACDM model with a "true" cosmological 
constant, the distribution of baryons can approximately be determined by the dark matter distribution, i.e., the 
observed power spectrum is proportional to 5 2 M . In dynamic dark-energy models as in the present case, however, the 
gravitational potential - which determines the baryonic distribution - is not only due to dark matter perturbations, 
but also due to the perturbations in the cosmological term. The total energy density perturbation 5 C is given by 

5 c^Pk±Pk=Pj± 5 c PX 5C (84) 
P P P 
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FIG. 1: Relative power spectrum d 85 1) as a function of k for Qmo = 0.3 (left) and Qmo = 0.8 (right). 
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FIG. 2: Relative power spectrum l|85|l as a function of k on large scales for !Jj/o = 0.3 (left) and Qmo — 0.8 (right). 



Only if the S x term is negligible, the observed spectrum can be related with 5^- The relative importance of the 
dark-energy fluctuations can be quantified by the expression 



PM 5 C M 



11-0 



3K 



MO Q 3/2 



MO 



d\nS 
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where we have applied equations (Q~9|) and ([66]) . In the integration of equation (|59|) wc have used scale-invariant, 
adiabatic initial conditions for k > at z = 1000, which is a very good approximation as discussed in subsection MI Cl 
above. In Figure [1] we show the ratio (|85[1 at present as a function of k for the range 0.01 Mpc _1 /i < k < 0.2 Mpc^ 1 h 
for two different values of iljuo- This ratio is much smaller than 1 in the entire region. Hence, we can see that 
perturbations in the vacuum term are actually negligible. The fact that (|85p decreases with k demonstrates that the 
dark energy component is perturbed significantly at the most at very large scales k ~ 0. But as Figure [2] reveals, the 
dark-energy perturbations remain smaller than the matter perturbations even on the largest scales. In Figure [3] we 
also show the ratio (|85[) &s a function of Omoj f° r two fixed scales. As to be expected, vacuum perturbations go to 
zero in the limits Omo = (de Sitter solution) and Omo = 1 (Einstein-de Sitter case). The result for larger k can be 
seen as a justification of the analysis performed in (l2| with px = 0. As already mentioned, a good fit of the observed 
spectrum was found in [12j for a present matter density parameter around 0.48. The corresponding best fit is shown 
in Figure [31 where the spectrum was normalized with the BBKS transfer function [13] for large k. 



It should be mentioned that the results of the present analysis are not in accordance with those presented by some 
of us in (l8j . where the vacuum term was perturbed and a fractional matter density near 1 was obtained. Some of the 
reasons for this discrepancy are the following, i) As initial condition, the ACDM-model based BBKS power spectrum 
at the redshift of last scattering has been used in [l8[ , something which is not completely justifiable since the present 
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FIG. 3: Relative power spectrum A85D as a function of Qmo f° r & — 0-01 (left) and k = 0.185 (right). 
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Log 10 k/h [Mpc- 1 ] 

FIG. 4: The best fit of the observed 2dFGRS power spectrum. The lines correspond to the BBKS transfer function (blue), the approximate 
ACDM best fit using the numerical method of [T3l (red) and the corresponding fit of the interacting model (violet). 

model is different from the ACDM model, ii) Only the conservation equations for the total fluid were considered in fl8| . 
but not the separate balance equations for the individual components. In this way, the non-adiabaticity of the system 
was not explicitly explored, iii) A separately conserved baryon component with a fixed relative density was included in 
[l8j . For this reason, the dark-matter density parameter could not vary in the entire interval [0,1], which led to some 
numerical difficulties. In the present gauge invariant treatment on the other hand, the non-adiabatic perturbations 
appear explicitly and scale-invariant initial conditions are implemented at very high rcdshifts. Our result corroborates 
the prevailing belief that in a fluid with cquation-of-state parameter uj = — 1 no significant perturbations should be 
expected on scales inside the horizon. 

As a final point we mention that for negligible dark energy perturbations px ~ the relative entropy perturbation 
Pnad in ([70)) with c 2 = — 1 reduces to 

Pnad ^ PX_ _ 1 1 - OmO ,ggv 

Pm P 2 1 - Q M0 + Q M0 a- 3 / 2 

Obviously, the entropic perturbations increase with the scale factor from a very small value for a C 1 to 1/2 in the 
long-time limit a ;§> 1. The present value is (1 — rij\fo)/2, leading to approximately 25% of non-adiabaticity. 



V. CONCLUSIONS 



We have investigated an interacting dark-energy model in which a vacuum term decays into dark matter linearly 
with the Hubble rate. The homogeneous and isotropic background evolution of this model is given by simple, explicitly 
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known functions of the scale factor. The two-component linear scalar perturbation dynamics was reduced to a single, 
second-order equation for the gauge-invariantly defined, comoving fractional density contrast of the dark matter 
component. Thanks to the specific coupling between the components, the dark-energy perturbation is determined 
by a combination of the dark-matter perturbation and its first time derivative. The presence of relative entropy 
perturbations makes the entire dynamics intrinsically non-adiabatic. However, as we demonstrated analytically, this 
non-adiabaticity is very small at early times. Consequently, for our numerical analysis adiabatic initial conditions 
could be used as a good approximation for k > at high redshifts. As a main result of this analysis we obtained 
that dark-energy fluctuations are negligible on scales that are relevant for cosmic structure formation. In other words, 
for our model dark energy does not cluster on small scales. This result provides an a posteriori justification of 
the analysis in (l2T | where dark-energy perturbations were neglected by assumption. This seems also to imply the 
conclusion obtained in [l2T | , that the matter-power spectrum is only consistent with a present matter-density parameter 
f^Mo ~ 0.48. This value is compatible (within 2a confidence level) with an updated joint background analysis of the 
A oc H model in [l4| • On scales larger than the horizon the contribution of dark-energy fluctuation increases but 
remains smaller than the dark- matter contribution even in the limit k — > 0. 
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